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1cm Abstract 

We calculate the radiative corrections to the decays r Mv^ and it lui, 
where the meson M is M = vr or i^" and the lepton / is / = e or /i. We perform 
a complete calculation, which includes internal bremsstrahlung and structure 
dependent radiation in the radiative decays and point meson, hadronic struc- 
ture dependent and short distance contributions in the virtual corrections. 
Our result for the radiative correction to the ratio r(r 7rz/T-(7))/r(7r 

/iz/^(7)) is 6Rr/^ = (O-ie^Q-JI) %. For the ratio T{t Kur{-f))/T{K 

IJ'i^/iil)), we obtain 6Rr/K = {^■^^^^ 26) "^"^^ completeness we have also 
calculated the ratio of the electronic and muonic decay modes of the pion. 

1 Introduction 



In the calculation of radiative corrections to semileptonic (semihadronic) decays such as 
r —>■ TiVj. and vr yuz/^, one faces three different problems [|T], 0]. As usual in radiative 
corrections, there are divergences, viz. first the infra-red (IR) divergences and second the 
ultra-violet (UV) divergences. The third problem, however, which is the central issue of 
this paper, is the treatment of the strong interaction. 

The IR divergences are removed as usual by considering either radiative decays with 
hard photons (eg. , r(r t^VtI) with > Eq) or inclusive rates for decays into final 
states with and without photon (eg. r(r — > ttz/t) + r(r — *• ttUt^)). 

The UV divergences are removed by renormalization. The decay rate of vr — fiu^ 
is proportional to the pion decay constant In principle is determined by the 
parameters of the standard model, but since we are not able to solve the nonperturbative 



regime of QCD, f^^ has to be considered as an additional free parameter which has to be 
extracted from experimental data and therefore has to be renormalized order by order in 
perturbation theory. On the other hand, ratios as 

_ r(r ^ nuril)) 
r(7r^/iz.,(7)) 

and 

_ r(7r ^ et/e(7)) 
r(7r /.^,(7)) 

are independent of and therefore can be predicted by theory. Technically the UV 
divergences cancel in these ratios. 

For a systematic treatment of the issue of strong interaction, a separation into different 
energy regimes should be made. If all momenta are very small compared to a typical 
hadronic resonance scale such as rup, the matrix elements are fixed by low energy theorems 
of QCD. In this low energy regime the pion behaves like a pointlike particle, and its 
interactions with the photon are determined by scalar QED. The high energy regime, on 
the other hand, is dominated by the short distance corrections to the weak interaction, 
ie. by photonic corrections acting at the quark level. These two regimes are not adjacent 
to each other, but rather there is an intermediate region which is dominated by non- 
perturbative strong interaction, viz. by the physics of hadron resonances such as the p 
and the ai particles. 

In a previous paper 0], we have calculated the corrections to Rt/w within a model 
with an effective point pion field. Defining the radiative correction 5Rr/-K by 

-Rt/tt = -Rr/lr(l + ^Rt/tt) (3) 

where 

^(0) ^ 1 ml {l-mllmlf 

denotes the prediction to the order 0(a°), we found a radiative correction 

<54^f ■) = +1.05%. (5) 

(P. M. denotes point meson). We added to this result another +0.17% arising from 
hadronic structure dependent radiation. So what is missing in are the hadronic struc- 
ture dependent effects in the virtual corrections and the short distance corrections. 

As for the short distance correction, it was shown in Ref. 0] that the leading 0(a) 
correction in the limit of a large Z boson mass, m\ — > cxd, is 



-Mo 



1 + ^ln^ 

TT /i 



Mo (6) 



where M.q denotes the Born amplitude, mz acts as a cut-off and /i is an unspecified mass 
scale characteristic of the process. The difficulty is of course to find identify the scale /i. 
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If the relevant scale is /ii for the decay r — ^ ni^ril) and fi2 for n /iz/^(7), the short 
distance contribution to the radiative correction of Rt/it is 

= ^ln^ (7) 

TT /ii 

and so only the difference in the scales fi2 and fii is relevant for the correction to Rr/n- 
In Ref. Marciano and Sirlin give an estimate for 6Rr/n which is based on the short 
distance contribution only, 

j^(s.d.) ^ ^(0) l + 2a/n\ia(mz/mr) .g. 

^/'^ 1 + f (a/vr) ln(m^/m^) + ^ia/n) \n{mz/mp) ^ ' 

is. d. denotes short distance) or 



2a . / ml/'^rn}J'^ ' 



5R[-t^ = - In '^^^^ ] = -0.98% (9) 



This result differs from the effective point meson result. But note that in the case of 
the pion decay this estimate extends the short distance physics to a very small scale of 
/i2 = m^/^my^ = 214 MeV. 

In a subsequent paper 0, Marciano included both the point meson and the short 
distance corrections to tt — li'i{'~f), matched at the scale rrip. For the r decay r — > 7rz/T-(7), 
still only the short distance corrections were included. Reexpressing his prediction in 
terms of SR^/n, we obtain from P, ^ 

5Rr/^ = (0.55 ±1.)% (10) 

where the ±1% is the author's estimate of the missing long distance corrections of 0{a) 
to T ^ 7rz/^(7). In a recent papper [Q, the authors have further improved the calculation 
of the radiative pion decay by including the leading hadronic structure dependent effects 
(both in the radiative decay and in the virtual corrections) and by including the leading 
two-loop effects. Using the pion decay constant /^r extracted from this calculation to 
predict the tau decay, they obtain a new prediction for the rate, which can be rewritten 
in terms of Rr/n as 

5Rr/n = {0.Q7±l.)% (11) 

where again the ±1% estimates the long distance corrections to r ^ nUri'y), which are 
still missing. 

Comparing the numbers in (|5|J5|, pi]| , pT] ) , it becomes clear that a complete and system- 
atic calculation of the full 0{a) corrections would be important in order to obtain a 
reliable prediction. This is what we will present in this paper. 

We have performed a systematic calculation of the radiative corrections to Rr/n 
which includes all relevant contributions. In the calculation of the loops, we separate 
the integration over the momentum of the virtual photon k into the long distance re- 
gion with < |/c^| ^fJ'lut into the short distance region with /icut^l^^l < ^'z- The 
matching scale ficut, which separates long and short wavelengths, should be of the or- 
der ficut ~ 0(1 GeV). The long distance part includes the effective point meson and 
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the hadronic structure dependent corrections. The latter is obtained by modifying the 
scalar QED coupling of the photon to the pion by vector meson dominance, and by 
adding loops which are proportional to the form factors Fy and which determine the 
hadronic structure dependent radiation. 

For the short distance corrections we first consider the leading logarithm, and second 
we give a complete calculation based on the parton model. 

Our paper is organized as follows: In Sec. 2 we briefly review the amplitudes for 
the radiative decays and the parametrizations of the two form factors appearing in the 
hadronic structure dependent radiation. In Sec. 3 we make some general remarks on our 
treatment of the virtual corrections, concerning the W boson propagator and the separa- 
tion into long and short distances. In Sec. 4 we calculate the point meson contribution. 
We have published the results for the radiative corrections in a model with a pointlike 
pion in 0] . Nevertheless we repeat them here for completeness and in order to give some 
more details and intermediate results, which are needed later in the calculation in order 
to combine the point meson correction with the other contributions. In Sec. 5 we consider 
the leading logarithm of the short distance contribution, and in Sec. 6 we calculate the 
non-leading corrections. In Sec. 7 we calculate the hadronic structure dependent loops . 
Then in Sec. 8 we explain how the different contributions combine to give the final result, 
which we evaluate numerically in Sec. 9. A summary and concluding remarks are given 
in Sec. 10. 



2 The Radiative Decays 

We consider the decays r Mi/t{'^) and M where the meson M is the pion 

M = TT or the kaon M = K, and the lepton / is the muon or the electron, / = /x or e. 
The Born amplitudes are given by 

Mo{r{s) ^ M{p)uM) = -GFVMfM[uuiq)]^l-Uris)] 



Mo{M{p) ^ l{sMq)) = GFVMfM[uu{q)l/lMs)] 



:i2) 



where 



V^ = cos 9c VK = sinec (13) 

The matrix elements of the radiative decays r — > Mv^-'-y ||^, |T0[ and M Iva ||Tl| 

can be written as the sums of two contributions, the internal bremsstrahlung (IB) and 
the structure dependent radiation (SD) 



M[t~{s) ^ Ur{q)M-{p)^{k)] = MlB + MsD 



(14) 



where 



MlB 
MsD 



-GpVMefMmr 



p ■ e 



+ 



s ■ e 



GpVMej 
V2 I 



u^{q)Yl-UT{s 



p ■ k 2s ■ k s ■ k 

F^^'\{k + pf 



t^kPp" 



rriM 
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rriM 



Similarly 
where 



M[M+{p)^l+{sMq)^{k)]=M'js + M'sn 



(15) 



(16) 



-GpVMefMm 



GpVMei. 
■{ le 



fiupcr 



\J-k- 2^k- 

F^^'\{k-pf\ 



e^kPp" 



+ 



Uu{q)l+{{p ■ k)i - (e ■ p)H)vi{s) 



Fr\ik-P? 
rriM 



(17) 



The IB part is fixed by the QCD low energy theorems. It is exactly identical to the 
one which would be obtained for an elementary pointlike pion field, with electromagnetic 
interactions determined by scalar QED. The IB amplitude does not contain any unknown 
parameters beyond the meson decay constant /m- The SD contribution, on the other 
hand, involves the two form factors Fy ^\t) and F^^\t), which describe the effects of 
non-perturbative strong interactions. Crossing symmetry implies that both the meson 
and the tau decays are described by the same analytical functions of the momentum 
transfer t, the difference being that in the r decay, t varies from m\j up to m^, and 







, m 



in radiative meson decay, t 
factors. 

For the pionic case M = vr, we gave 



In Ref. [|T0| we have parametrized these form 



FW(t) = FA-^(0)BW,,(t) 



and 



Fa'' (0) has been measured in radiative pion decay |T2| 



Fl^\t) = 4"^(0) [BWp(t) + (TBWp.(t) + pBWp.(t)] - 



a + p 



fP{o) 



-0.0116 ±0.0016. 



whereas Fy\o) is related to the axial anomaly and predicted to be 



Fr(o) 



+0.0270 



(18) 



(19) 



(20) 



(21) 



Note that here we use a relative sign s := sign{f.,rFv{0)) which is positive in our conven- 



tions, s = +1. In fact in 10 1, we used s = — 1. However, as we have explained in [13 



we 



now beheve that s = +1 is the physical choice, which we will therefore use throughout 
this paper. 
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The Breit-Wigner resonance factors BWx(t) are normalized according to BWx(O) = 1 
and involve energy dependent widths 



or constants widths 



BWx(t) 



BWx(t) 



X 



— inixTx 



m\ — t — imx^ 



(22) 



(23) 



X 



For the ai resonance, an energy dependent width based on the three pion and the pion-rho 
phase space has been calculated in WM, 



ra,(t) 



/ 2 \ ''l 



(24) 



with 

9it)- 





A.l{t - 9ml)^[l - 3.3(t - 9ml) + 5.8{t - 9ml)^] 
t(1.623 + 10.38/t2 - 9.32/t^ + 0.65/t^) 



if t < 9ml 

if 9ml ^ ^ < i^p + i^-n 
else 



(all numbers in appropriate powers of GeV) 

For the p, p' and p" resonances, an energy dependent width 



(25) 



rx(t) 





mx_ 



t — Ami 



m\ — Amt 



X 



if t < Ami 
else 



(26) 



can be derived from the P-wave two body phase space. 

We fixed the two parameters a and p in by employing four constraints, viz. the 
QCD theorem on hm FT,^{t) [|15| (where Fy\t) = —mT,F^j{t)/^/2 [0), the measurement 
of the slope F'(t = 0) and the widths Fp^^^ and Tpi^T,^. The result was 



a 
P 



+0.136 
-0.051 



(27) 



We also compared to a dipole parametrization (a = 0.0584, p = 0) and to a monopole 
form {a = p = 0). 

For the kaonic form factors (M = K), consider the parametrizations 



V 



it) 
(t) 



4'^)(0) [BW K^t) + CTKBWK^it) + pKBWK*"{t)] 
Ff\0)BWK,il270){t) 



1 



1 + aK + Pk 



(28) 



(where K* = K*{892), K*' = K^(1410) and K*" = is:*(1680)). In |Tg we used the 
monopole parametrization ax = = of the vector form factor and constant widths 
for both the K* and the Ki. 
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Flavour symmetry implies the following relations for the form factors at t = 0: 

fW(0) = !^F^)(o) = +0.0410 

4^)(0) = !!^4-)(0) = +0.0955 (29) 



The measurement of K i^il, on the other hand |12], gives 



F^^^ + Pf^ = +0.148 ±0.010 

F^^'^-Ff^ e [-0.3,2.2] (30) 

Taking the vector form factor from flavour symmetry and the anomaly and calculating 
the axial one from the measured sum and the value for Fy^\o), this results in 

Ff\0) = —M'^VO) = +0.0525 ± 0.010 

4'')(0) = —4'')(0) = +0.0955 (31) 
These are the values we will actually use. 

3 General Considerations on the Treatment of the 
Virtual Corrections 

The momentum dependence of the W propagator 

^ ""'^ (32) 



is determined by the familiar Feynman cut-off function m^/(m^ — g^). Thus we can 
as well use a local interaction with an UV cut-off equal to mw in the calculation of the 
virtual corrections. In fact, it has been shown by Sirlin [ITS that the radiative corrections 



of the order Gpcn calculated within the full standard model with a single Higgs doublet are 
equal to the photonic corrections calculated with the local V — A interaction and a cut-off 
equal to rriz, except for very small contributions of the order asGpCt- While the photonic 
corrections are identical to those computed in the local V — A theory with an effective 
cut-off equal to m^y, the non-photonic corrections lead to an additional contribution 
which depends on \n{mz/Tn\Y) and 6^. In the simplest electroweak model with a single 
Higgs doublet, where cos 9w = rriw/mz, the photonic and the non-photonic corrections 
combine to give a result which is identical to the photonic correction obtained in the local 
theory with the cut-off equal to mz- As we will show, the actual value of the cut-off does 
not matter for our results on Rr/n, because it cancels in this ratio. 

Thus we will use a local V — A interaction and dimensional regularization, ie. the loop 
integrals are evaluated in 4 — e rather than in 4 space-time dimensions. We can translate 
our results into a momentum space cut-off by the replacement 

A - In — In — I 33 
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where 

2 

A := - -7s„;er + ln47r (34) 

The virtual corrections fall into three classes, see Fig. ||. The first class includes only 
a single diagram in which all couplings are known, and the third class is identical for 
r TCUr and vr fiu^, so these diagrams drop out in the ratio Rr/n- Therefore the most 
important class for the calculation of Rt/-k is the second one. 

Now we want to separate the integration over the momentum of the virtual photon 
into two regions with small and large k'^, respectively [T7|. We achieve this by splitting 
the photon propagator: 

+ 13. 7T- (35) 



"long distance" "short distance" 

Obviously the first part is important only for \k'^\^fJ'lut the second part only for 
^fJilutt so indeed they correspond to long and short distances, respectively. And so 
the photon propagator is divided into two parts, a regulated photon propagator with an 
effective cut-off ^cut and a massive photon propagator with mass ^icut- The scale ^cut 
separates long and short wavelengths and should be of the order ^^ut ~ 0(1 GeV). 

A photon with short wavelength {\k'^\^ l^'iut) resolves the quarks in the pion. It in- 
teracts with the quarks which come into being in the initial process r v^-ud, which 
subsesquently hadronize to form the pion (see Fig. |^). A photon with long wavelength 
(l^^l ^A'-cMt)' other hand, has a small resolution and interacts with the pion as a 

whole or perhaps with some hadronic resonances (see Fig. ^. So according to Eqn. ( pS]) 
we have to calculate the short distance diagrams such as in Fig. § with a massive photon 
propagator and the long distance diagrams such as in Fig. |^ with a regulated photon 
propagator. 

In the next section we will start with the calculation of the corrections in a point pion 
model. This will give a first estimate, and it is part of the complete calculation. 



4 Point Meson Contribution 

We will now calculate the corrections to r — tti/t- and vr — > /xz/^ in a model with an 
effective point pion field, using a generalized photon propagator with a mass 



P — k"^ — im? 

(A denotes a small IR regulator mass, which in the end of the calculation is put equal to 
zero, whereas m denotes any finite mass.) It will become clear below why the considera- 
tion of such a generalized photon propagator is useful, see Eqns. (|TTD|)-(|TT3]). 

For the case of the pion decay (and for = A^), these corrections have been calcu- 
lated in 0. Our calculation, however, differs in some technical details from that in 0. 
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First we use dimensional rather that cut-off regularization. Second Kinoshita replaces 
the vector-minus-axial vector current interaction 

GpcosdcWini - 75)^.](^5^ - eA^)$. (37) 
by the scalar-minus-pseudoscalar current interaction 

G^cos^cAmi°)[^/(l -75)^.]$. (38) 
where mf'^ denotes the bare lepton mass. We have performed the calculation both using 



Eqn. (pTf) and using Eqn. (|3^) and we have checked that they give identical results. Here 
we present the results using the V — A form of Eqn. (p7|). 

There are six Feynman diagrams for the virtual corrections to r ^ vrz/^ , see Fig. ^, 
and of course there are six similar diagrams for vr — > fiu^. The last diagram SAix actually 
vanishes after pion mass renormalization. 

The ratios of the virtual correction amplitudes over the Born amplitudes for both the 
tau and the meson decay can be expressed by the same functions of mf and m^j, if we 
use a general lepton mass mi which denotes rrir, rrifj, or me, respectively. The result is 

^(m^mL,m^) = ^|2< - < + 2mfCo - 2mLCi| 
^(m^mL,m^) = ^jl - 45^ - 2i?l| 

^(m^mL,m^) = ^jl + i^j + [4< + 2i?i']| (39) 

where the B\, and Ci are given in terms of the standard n-point functions B and C: 

Ci = Ci{ml[,m'f,0,mM,m,mi) 
B\ = Bi{m^,mi,m) 
Bf = B,{ml,mM,m) (40) 

(see App. ^ for our conventions). Thus 

SMpML. 2 2 2n ^^SMi 

+m] [45;' + 2B[) + 2m^Co - 2ml,Ci \ (41) 
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("PML" denotes "point meson loops".) This is the form which will be used in the 
calculation of the complete radiative correction in Sec. |^. 

To obtain the radiative correction within a model with a pointlike pion, we now 
assume = and calculate the n-point functions. The result is 



5 M PML I 2 2 x2n « J 3 3 m^,^ 



^Inn + 1 



r 



rriM 



1 — r, \ 1 — r. 



+ 2- -L(lnrz)2+ (5-^^)lnn^ (42) 



where 

r, := ^ (43) 

rriM 

and is the mass scale of dimensional regularization. 

By interference with the Born amplitude, the virtual correction of the decay rate is 
of course 

To ~ Mo ^^^^ 

To these virtual corrections the integrated decay rate for internal bremsstrahlung must 
be added. We divide this rate into "soft" {E^ < Eq) and "hard" {E^ > Eq) corrections, 

where Eq is assumed to be small, ie. Xq := (mT-/2)_E'o ^ 1- For the r decay the results 
are 



2 



+ 4i±4^ In TM ln(l - rl,) + 0(xo) [ (46) 



and 



+ (1 - 1-^) - 4 ln(l - rl,) + 2^-±^U,{rl,) 



M ^ ' M 



where 



+ 0{xo) (47) 



tm ■= (4J 



10 



For the meson decay the sum of soft and hard bremsstrahlung is 



a f, /1 + r? \, A l + rf 



4 Inn + ij ln(l - rf) + J^,^, Inr, 

.1 + 2N 27 3 1 

4^-iLi2(l-rr ' ^ 



1-rf •«/ ■ 4 21 -rf 



In these formulae, Li2 denotes the dilogarithmic function 

ln(l - t) 



Li^ix) = - /' 

^0 



dt 



t 



Now we add up virtual, soft and hard photonic corrections as obtained in 
meson model. Writing the radiatively corrected rates as 

r (channel) SF 



ro (channel) 



1 + channel^ 



we obtain 



(m ^ /i.,(7)) = - A + In + 61nn + - - -TT^ + /(n) 



'^'V 27r[ 2 2 '43 

and 

po V ---^v'V 27rl 2" ' 2"^ //2 ' 4 3 



<^r/ . ^\ tt f 3^ 3, m2 17 2 2 , ^" 



where 

fin) 



4([^lnn-l)ln(l-rf)-^li^lnn 

+4l±4Li.(rf)-^f^ + i.^ 
1 -r,2 1 -r? V2 3 



and 



givM) = 4 (\^lnrM - l] H^-rl) - ^ InrM 



l + '^ifT- / 2 A , /3 4 
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We have written the corrections in such way that /(O) = g{0) = 0, so from Eqn. ( pSj ) 
we get the well-known lepton mass singularity of the total radiative correction to pion 
decay as 3a/7rln(m;), whereas in the radiative correction to the tau decay the meson 
mass singularities cancel according to Eqn. (0). Our results are a nice application 
of the Kinoshita-Lee-Naunenberg theorem, which states that mass singularities cancel 



in inclusive decay rates g, |T9|. In the case of the pion decay the Born amplitude is 
proportional to mi and therefore a logarithm In mi is allowed in the radiative correction, 
while for the tau decay the Born amplitude is not proportional to m^, and therefore a 
logarithm In m^^ is forbidden in the radiative correction. 

The ultra-violet divergences of the corrections to the tau and the meson decays are 
equal and they cancel in the ratio Rt/m defined by 

with the finite radiative correction 



V ' J PM 27r I 2 m%j mu 2 \mr J \mM J ] 



So within this model with an effective pointlike meson ("PM"), we end up with the 
results H 



{5Rr/.) 
{dRr/K) 



= +1.05% 

PM 



PM 



+ 1.67% (59) 



Note that this result differs from the Marciano-Sirlin estimate quoted in the introduction. 



6R'',f = -0.98% 



t/tv 
iM.S. 







6R^'^-^- = -0.53% (60) 

This is mainly due to the fact that the point meson and the short distance corrections 
do not have the same UV divergences. 

From Eqn. (^) the radiative correction to the ratio of the electronic and the muonic 
decay modes of the pion 

r(7r ez/e(7)) _ml f m^ 
r(7r /iZ/^(7)) ml \m: 

can also be calculated in the same way and we obtain 



Re/, = T^ T^ = ^ ^ 1 \^+^Re/,] (61) 




{5Re/i) 



PM 




-3.93% (62) 



a result derived long ago by Kinoshita 
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5 The Leading Short Distance Logarithm 



In Ref. 1^ the author shows that for any semileptonic weak process, the leading 0{a) 
correction in the hmit of a large Z boson mass, m| — oo, is 



^0 



1 + ^ln^ 



Mr 



(63) 



where A4o denotes the Born amplitude and yU is an unspecified mass scale characteristic 
of the process. 

We will now rederive this result and then show that we are able to fix the scale /i in 
terms of the lepton mass mi and the scale ficut- Consider the short distance corrections 
to the amplitude for the initial weak process r Urud 



.Gf cos 6', 



c , 



V2 



(64) 



The six Feynman diagrams which give the radiative corrections to this short distance 
amplitude are shown in Fig. ^ In order to find the term of the order 0{a\nmz), we 
have to calculate the UV divergence of the short distance correction 



6 A = 6Aa + 6Ab + ... + 6A 



f 



(65) 



The calculation can be simplified by using the Landau gauge In this gauge the 

amplitudes corresponding to external line renormalization, 6Aa ■ ■ ■ SAc, are UV finite, as 
are SAe and SAf, where the photon loop connects an ingoing with an outgoing fermion 
line. In the Landau gauge, the only UV divergent amplitude is 6Ad- 



SA, 



(2n 



(k 



2^3 



dap 



ka kj3 



—AAo + ... 
Zn 



+ ... 

(66) 



where the dots indicate terms which are UV finite. 

Now replacing A by the UV cut-off lnm|, we obtain 

5A=-\ti — Ao + ... 

71 fl 



(67) 



where fi is some unspecified characteristic mass scale which must be introduced by di- 
mensional arguments. This is Sirlin's result. 

If we replace the photon propagator by 1/(A;^ — yUcut) if we do not neglect 
the lepton mass, we obtain 



6A 



a 



1 



In ^cut In 



A^ 



(68) 



We will now apply this result to the radiative correction to -Rt/tt- In addition to the 
short distance correction which contains the above logarithm, there are the long distance 
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corrections associated with the integration over = . . ■ fJ^lut^ which have to be added. 
Of course these long distance corrections have to be computed using the effective pointhke 
pion (and hadronic resonances) and give rise to other logarithms which in fact turn out 
to be more important in Rr/n that the short distance ones. 

So we will now combine the long and the short distance corrections. First we have 
to discuss which scale is a good choice for ficut- The effective point meson is a good 
approximation only if all momentum transfers squared are small compared to m^, and 
so for the long distance part one should require <^ m^. The short distance part, 
on the other hand, uses asymptotically free quarks which one would believe in only for 
/^c«t > (1 . . .2 GeV)^. A compromise between these non-overlapping regions would be 

2 2 

The long distance corrections are to be integrated over fc^ = . . . /i^^^ , which is taken 
into account by using an UV cut-off ficut = for the point meson results of Sec. H (ie. 
by the replacement A ln(m^//i^) in Eqns. (^) and (|5^). Thus we obtain 



short dist. loi^g dist. 



, / n\ 2a;, mz a ( 3 . ml ^, 11 2 „, 



and 



2.22% -1.02% 



short dist. 



— (r^7rz/^ 7 l = ^ m^ln m In 

1 ^ ^ ^ vr — y m,- nip 



1.74% 

long dist. 



a (3, ml 17 2 ^ , ,1 



0.03% 

This modifies the point meson results in the following way: 

e-T-. /e-T-. \ 2a ml , 

V ' / PM vr mi — mi m^ 



Similarly 



-0.48% 

0.57% (71) 



/rr-. \ 2a; ml . mx* 

5Rr/K = [^Rrl.)^^, + — In- 



P^i TT mi — mx* mj 



-0.43% 

1.24% (72) 
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The prediction for 8Rr/-,r obtained in this way by matching the leading short distance 
logarithms with the effective point pion correction, SR-t/tt = 0.57%, differs strongly from 
the first Marciano-Sirlin estimate ||^, SRr/n = —0.98%, which is based on the short 
distance logarithms only. Note, however, that it is numerically close to the predictions 



in H, see Eqns. dlQilTlD above. The estimate in g (and up to very small additional 
corrections, the value in as well) includes long and short distance corrections as in 
Eqn. dl) for the pion decay. For the tau decay, the authors only include the short 
distance correction, estimated using Eqn. (|63| ) with fi = 171^, ie. they estimate the distance 
correction to be (2a/7r) ln(m^/mT-) = 1.83%, which has to be compared to the 1.74% 
calculated in Eqn. {^^)■ 

And so, what is really missing in [^, §], as compared to our estimate in Eqn. (|7l|), are 
the long distance corrections to the tau decay. In || the authors only estimated their 
possible size, resulting in a ±1% uncertainty. As we have calculated, the long distance 
corrections in the model with an effective point meson happen to be extremely small 
(0.03%), which explains why our estimate in Eqn. (^) and the estimates in |] are 
quite similar numerically. 

However, these estimates still cannot be considered as safe for the following reasons: 

1. The value ficut = is too large for the point meson approximation in the long 
distance part and too small for the assumption of almost free quarks in the short 
distance part. Therefore the range of validity in the long distance part should 
be extended to 1 ... 2 GeV by including vector meson resonances and then Hcut = 
1 ... 2 GeV will be a good value to use. Indeed we will show below that the vector 
meson effects in the loops change the final result considerably. 

2. For the long distance corrections to the tau decay, the UV cut-off ficut is smaller 
than rriT-. Therefore terms proportional to ml/fil^f are missing in Eqn. (|70|). 



3. The large scale mz appearing in the short distance logarithm cancels in the ratio 
Rt/tt, leaving logarithms of comparable scales (Inm^ and Inm,-). Therefore it is not 
obvious that non-logarithmic contributions in the short distance correction can be 
neglected. 

So in the following sections we will improve the calculation of the long distance cor- 
rections by properly using the regulated photon propagator of Eqn. (|35|) and by including 
vector meson resonances, and that of of the short distance corrections by going beyond 
the leading logarithm. 



6 Short Distance Beyond Leading Logarithm 

In this section we will present a complete calculation of the short distance corrections 
beyond the logarithm which is leading in the limit m| 00. This leading logarithmic 
contribution to the correction 6 A of the short distance amplitude Ao{t — >■ VT-ud) could 
be written as a factor times the Born amplitude, 5 A = C x Ao + . . .. Then the correction 
6Ai to the exclusive rate A^o(t ^ ttz/,-) involves the same logarithm, 6Ai = C x Aio + ■ ■ ■■ 
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The situation, however, is more involved for the complete result for SA. Firstly, in general 
it can not be written as a number times but rather it involves other operators, and 
secondly 6 A will depend on the relative momentum of the two quarks. The first problem 
can be solved by projection on the = 0~ state of the two quarks, and the second 
problem by use of the parton model. 

Let the short distance amplitude for the decay 



(73) 



of the tau into neutrino and an quark- ant iquark pair with relative momentum 21 be given 
by A{1) (and similarly for the decay into us, for the time being we consider the case of 
non-strange decay into the pion only). In a frame, where the pion is moving with infinite 
momentum, / is proportional to p, 

II 

(74) 



/ = —p 



and A = A{u) . Then the amplitude of the exclusive tau decay r 

M{r- ^n-Ur) =T{A) 

where 



TTi/^ is (see Fig. ^) 
(75) 



TiA) 



-1 



rrir, 



A{u)ui, 



1 



u 



P IbUd 



l+u 



-p 



(76) 



Here $7r(M) = $7r(— w) is a symmetric parton distribution function (The numerical factor 
is just conventional. Our treatment follows closely that of Refs. |]21|, 0). The sum ^ is 
over the spins of the quarks. For the quark mass we assume isospin (and in the case 
of the kaon, S't/(3)-flavour) symmetry. 



nir. 



(77) 



The limit — is implied. Note that all relevant operators in A involve an odd number 
of Dirac matrices between the Ud and the Uu Dirac spinors, so the sum is proportional to 
iTLq, and this rUq cancel the rriq in the denominator. 



In the case of the Born amplitude M.^, Aq is 



Ao = -i 



Gp cos 6c 



and so 



■-Gp cosOc / du 
8 j-i 



T{Ao) 
+1 $.(«) 



)tr 



1-u 



TTlr 



75 



l + u 



--Gp cos Or 
2 



+1 



1 



U^p/'J-Ur 



(79) 
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where the last hne follows from the definition of /tt- So we have rederived the sum rule 

-1 



|23 



du(f)^{u) = 



(80) 



After these preparations, we can now calculate the full short distance corrections 
5Ai for i = a . . . f given by the Feynman diagrams in Fig. ^, using a massive photon 
propagator with — > ~ lAut)^ Feynman gauge, and neglecting the quark masses. 
According to Eqn. (fTSD, the corrections 5Ai induce corrections 5A^j(r Mv^) to the 
rate for the tau decay 



5Mi{T Mur) = T{6Ai 



{i = a...f) 



(81) 



(and analogously for the pion decay). The results depend on the unknown distribution 
function ^m{u) and are given in the form of integrals over u 



^■^if 2 2 2 ^ 



+1 



cut) 



with weight functions Wi{u). Define the total short distance correction amplitude 



^Msd, 2 2 2 \ 

and the sum of the weight functions 



ESAii 2 2 2 \ 
— — (m;,mA^,/i,J 
i=a -^0 



(82) 



(83) 



w{u, m], ml^, = ^ Wi{u, , m|^, /i^^J 



(84) 



Then the short distance correction {SRr/M)sd is given by 



{SRr/M)sd — 2 



6M 



sd/ 2 2 2 ^1 

[m^, m^, ^cut) 



6M 



Mo 



sd/ 2 2 2 \ 
'^/i' f'cut) 



— / du(j)M{u)[w{u,m^,m^^j,fi^^^) -w{u,m m^,n^^^)] 



(85) 



Note that Wa, Wb and We are identical for tau and pion decay, ie. they do not depend on 
the lepton mass. Using furthermore the symmetry of the distribution function (pM under 



u —u, we can write 



[SRr/Al) sd 



2/m J-i 



du(f)M{u)rr/M{u) 



with 



w{u, ml, m^, fil^t) - w{u, ml, ml^, /i^^J 
where w is the symmetrized sum of the weightfunctions Wc, Wd and wj: 

\ [wc{.u, mf,mlj, + Wd{u, mf,ml^, fi^ 

+Wf{u, mf, mlj, fil^t) - Wc{-u, mf, ml^, 
-Wd{-u, mf, m\^, fil^t) - Wf{-u, mf, mlj, fif^^) 



(86) 
(87) 



w{u,mlml^,fil^^] 



(88) 



17 



Similarly the short distance correction to Re/fi is 

3 '•I 



me/,) 



sd 



2fM J-i 



with 



re/^,{u) := 2 
The relevant weight functions are given by 



w{u, ml, ml, - w{u, ml, ml, ^i^^^] 



a 



— n - 2fio - 2fii - 4mf[5; - B'q] 



57r ^ 

Wd{u,ml,m]^,^il^^) 
a 

+ [(1 + ufmij + (1 + u)mi]Cf + [(1 + u)mij + (3 + u)mi]C^ 
+ {l + u)mlC^^-2\ 



L f ivi I I cat/ 

;|l6Co^o + (1 + uf^liCii + ^rnlC^^ + 2(1 + n)(m^ + ml)C- 



(89) 
(90) 



where 



a 

12^ 



|4Ci(o + (1 - ufmlfil, + 2(mi, + m^C^ + (1 - w)(3mi, + m^C^ 
+ [(1 - ufmli + (1 - n)mf]C/ + [(1 - u)mli + (3 - u)ml]Ci 
+ {l-u)mlcd (91) 



B['}{ml,^cuum) 

( {I + Uf 2 2 

C I ^ mlj, — z — mf 



C 



4 

\l-uf 



2 

1 +M 



1-n^ 
4 

1-^2 



mM,ml,ficut,0,mi 



-mM,mi,Hcut,0,mi 



(92) 



In Figs. ^ - 1^ we have plotted the symmetric functions rT-/7r(M), r^iK{u) and re/^{u) 
for three different values of /icut- It turns out that these functions change only very little 
while varying u from u = 0, where the distribution function (p^^ is presumably peaked, to 
u = 1. Therefore we can approximate the weight functions by their values at m = 0: 



To 



^mf, m\j, /icut) 



3 

-— ti;i(0,m;2,m^,/x2^J / du<j)M{u) 

ZJm "'-1 



Wi{0,ml,mlj,fil^t) 



(93) 
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and then 



3 /"i 

{5Rr/M)sd = TT^^r/MiO) / du(j)M{u) 



= r^MO) (94) 
It is clear from Figs. - ^ that this approximation induces an uncertainty which is well 

below tool% for ^-/-' well below toof% ^r/K and well below ^qqqq% for Re/^,■ 
Note that we can make these statements without any specific assumptions on the form 
of the distribution functions 0,r and (px- 

In Fig. ^ we compare {SRr/n)sd as a function of ficut as given by Eqn. ( pi]) with the 
estimate from the leading logarithms 

[oRr/njsd ~ 2~ 1^ (95) 

IT ml- ni^t mr 



We find that Eqn. (^) gives an excellent approximation to Eqn. ( p4D and that the non- 
leading contributions to (5-Rr/7r)sd are very small. 



7 Hadronic Structure Dependent Loops 

For the real photon emission we have calculated the hadronic structure dependent ef- 
fects in [|10|, cf. Sec. @. In the long distance virtual corrections there are two different 
effects associated with hadronic structure. On the one hand, the photon emitted in the 
radiative decay by hadronic structure dependent radiation (SD) as a real photon could 
be reabsorbed either by the lepton or by the pion (kaon), see Fig. |11|. If in the radiative 
decay (fc^ = 0) there is a SD amplitude proportional to Fy and Fa , the corresponding 
hadronic "structure dependent loops" (SDL) must also be there for sufficiently small k'^ 
(long distance). Actually for the last two diagrams, where the photon couples to the me- 
son, the respective corrections to tau and to meson decay will be identical and therefore 
they will cancel in -Rt/tt- So for simplicity we will not consider them. 

On the other hand, according to the notion of vector meson dominance, the photon 
does not couple directly to the pion (kaon) but rather through a p intermediate state, 
see Fig. |12|. Similarly in the diagrams proportional to Fy and F^ the photon might 
couple through vector meson dominance as indicated in Fig. In the remaining of 
this article these modifications of the photon couplings will only be called "vector meson 
dominance" of the respective diagrams, either of the point meson loops (PML) or of the 
hadronic structure dependent loops (SDL). The name "hadronic structure dependent" 
will be used only for amplitudes such as and 5A^7 which are proportional to the 

form factors Fy and Fa- This is of course just a naming convention; the vector meson 
dominance of the photon coupling is also an effect of hadronic structure. 

We will discuss the modifications due to vector meson dominance of the photon cou- 
plings later and start with the calculation of SAig and SAi^ using the generalized photon 
propagator with mass 
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and no vector meson dominance in the coupling of the photon to the vector meson and 
the pion (see Figs. [Tl| (a) and (b)). 
The result for 6J^e is 

-j^{mi,mM,My,m ) = — J^mMfM ^ ^ - "^m)(<^22 + ^12)) (97) 

where 

Cij = Cij{ml^, 0, m^, m, My, mi) (98) 
My denotes a - possibly complex — vector meson mass 

= - imvTv (99) 

A comment on the treatment of the meson propagators for the p, p' , ai, ... particles is in 
order. In the calculation of the radiative tau decay we used sophisticated Breit-Wigner 
resonance factors with energy dependent widths. However, the difference between fixed 
and energy dependent widths is very small after integration over the spectrum. In the 
virtual corrections one has to integrate over all possible loop momenta anyway. For the 
virtual corrections, we will therefore only use Breit-Wigner resonance factors with fixed 
widths 

^^^^^ ) - - ;;tv 

or even with real masses My = rriy. 

Similarly the amplitude SAi^ is equal to 

Mo AtT y/2mMfM ^ 

+ (2m2 + 4m^)Ci2 + (m^ - m^,)Ci - l} (101) 

with 

C... = C..Xml„0,mlm,MA,mi) (102) 

In the case of the meson being the pion, our standard parameterization of the form factor 
Fy\t) in the radiative decay was a tripole dominated by resonances p, p' and p" with 
relative strengths a and p (see Sec. The same parameterization will be used here for 
the structure dependent loop 6Me , resulting in 

^Mv f 2 2 2n 1 fSMe 2 2 ,^2 2n 

Mo 1 + a + p^ Mo ^ 

6M 



M 







Mo 

^(m^m^,m^) = ^(m^ Mf^, m^) (103) 
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and for the kaon we define 



Mo I + aK + Pk ^ Mo 

6MIq , 2 2 71 /r2 2\ 
. 2 2 T\ r2 2\\ 

+ PK -r-j-{mi,mj^,Mj^.„,m )j 



Mo 

^-^{mlml^m-) = ^K^mi,Mi,m^) (104) 



Here again we use the convention that capital letters for masses indicate possibly 
complex masses, 

= "^ai - imaj^ai = - iniKj^ (105) 

and so on. In terms of the rate the corrections are then given by 

^^SDL, 2 2 2\ oT> {^-^V, 2 2 2n , ^Ma/ 2 2 2\\ l^(\a\ 

(m,,m^,m ) = 2Re | (m^ , m^v^ , m ) + -^-(m;,mA^,m )\ (106) 

Now we will discuss diagram by diagram whether the coupling to the photon should be 
modified by vector meson dominance. Consider the effective point meson graphs in Fig. ^. 
In the coupling 77r7r is determined by the electromagnetic form factor of the pion, 

which is well know to be dominated by the p vector meson (see Fig. |T2|). The diagrams 



bM-i and bMi^ cancel in the ratio -Rt/tt? so we will not include vector meson dominance 
here. In diagram bM^ the photon couples to the tau only, so the only diagram which 
remains to be considered apart of biM\ is the diagram bM2 with the seagull coupling of 
the photon to the weak interaction vertex. Here it not clear at all whether this graph 



should be multiplied by a rho Breit Wigner (compare Fig. |T4D. However, care must be 
taken to insure gauge invariance. The individual diagrams are not gauge invariant, but 
their sum is. So the modification of bM2 must be made in such a way that the sum of 
the diagrams is gauge invariant. This imposition of gauge invariance determines that the 
correct modification of the diagram bM2 is given neither by the multiplication with 1 
(ie. no vector meson dominance, VMD) nor by the multiplication with BWp (complete 
VMD), but by the multiplication with 

2BWp(fc2) - 1 (107) 

We will call this below the "seagull type VMD" . 

Next consider the corrections bM% and bM-j. In these diagrams the coupling of the 
photon could be dominated by the uj and the p mesons, respectively (see Fig. |l^). We 
do not have any experimental information on whether or not these couplings are actually 
dominated by vector mesons. (An experimental test could be made by measuring the 
e^e~ invariant mass spectra in the decays p — * 'ne^e~ and a\ TTe'^e~ , respectively.) 
By extrapolation from the experience with other hadronic couplings of the photon one 
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could expect vector meson dominance here as well, but in order to be unprejudiced we 
will below consider both possibilities, complete VMD and no VMD. 

Now according to Sec. ^, in the long distance diagrams a regularized photon propa- 
gator 

is to be used. If the photon additionally couples via vector meson dominance with a 
vector resonance mass M^, the simple photon propagator has to be replaced by 



fc2 - A2 A;2 - A2 -k^Ml- 
Assume a long distance diagram has been calculated with the photon propagator 

1 



^2 _ ^2 



(110) 



and without vector meson dominance coupling of the photon, with the result G{m^). 
Then the following replacements have to be made in order to get the correct answer: 

• long distance diagram without VMD 

G(m^) - G(A^) - (111) 

• long distance diagram with usual VMD 



G(m^) - G{\') - Gif^lJ + [Gifi'tJ - G{Ml)\ (112) 

long distance diagram with seagull type VMD 



2^2 

l-^cut — Mj^ 



G(m2) ^ GiV) - G(/iiJ + .2 Gif^Lt) - G(M2) (113) 



8 Complete Radiative Correction 

We write the complete radiative correction as 

SRr/M = {^Rt/m)cPM+ {5Rt/m)vMD(PML)+ {5Rt/m)hSD 

+ (SRt/m) VMD(HSD) + i5Rr/M)sd (114) 

using the naming conventions 

• "CPM" = cut point meson, ie. the long distance correction due to an effective point 
meson, including real and virtual photonic correction, where the loops have been 
calculated with the regulated photon propagator with the cut-off scale fi^^^ 
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• "VMD(PML)" = vector meson dominance of the point meson loops, ie. the correc- 
tion of the CPM result due to the vector meson dominance of the meson electro- 
magnetic form factor, 

• "HSD" = hadronic structure dependent, ie. the correction due to the diagrams 
proportional to Fy and F^, including real and virtual corrections, 

• "VMD(SDL)" = vector meson dominance of the structure dependent loops, ie. the 
modification of the HSD result due to vector meson dominance coupling of the 
photon in these diagrams, and finally 

• "sd" , the short distance correction. 

For the integrated rates of the real photon emission we divide into pure internal bremsstrahlung 
(IB) and the rest (SD -|- INT), viz. the sum of pure structure dependent radiation and 
of the interference between the internal bremsstrahlung and the structure dependent 
radiation. 



And so 

{5Rr/M)cPM 



ro(r ^ Mv^) To{M ^ fiu^) 

— (m^,mJ^^,A ) — (m^, m^^, 

i i 

STpML, 2 2 \2\ , STpML, 2 2 2 



and 

{SRt/m)hsd 



i i 

TsD+iNT{r Mt/^7) _ rsD+iNT{M iiu^-f) 
Voir ^ Mur) ro(M ^ iiu^) 

, ST SDL, 2 2 \2\ ST SDL, 2 2 2 \ 
i i 

-^^^{rnl,ml^,X'^) + ^^^{rnl,rnl^,^l^t) (116) 
i i 

For the point meson loops the vector meson dominance is taken into account by 

/r> \ ^VMD(PML) ,0 2 \ T vMD(PML) , 2 2 \ /n^x 

{Rt/m) vmd(pml) = ^("^T> t^m) F ^(^M' ^m) (117) 

i i 

where according to the last section 

i 1 CTi pi I 



t4ut - Mj 



P L 



23 



5 Ml 



l^cut 



Mo 
5 Ml 



{mlml,Ml)-2^-^{mlml,M^ 



Mo 



Mo 



l^cut 



+^1,2 



P' L 

5 Ml 
5 Ml 



Mo 



mlmlj,Ml)-2 



5M 



Mo 



Mo ^■-"•-^-^'--^^ Mo 



(118) 



where the parameters ai and pi describe the relative contribution of the p and the p" 
in the electromagnetic form factor of the meson (compare a and p in the vector form 
factor Fv). Note that for the case of the kaon, ruM — ttir , this assumes complete U{3) 
flavour symmetry in the vector meson sector, Mp = M^^ = M^. Otherwise the relative 
contributions of the p", the cv and the $ in the electromagnetic form factor of the kaon 
would have to be considered. 

The vector meson dominance in the "hadronic structure dependent loops" (SDL) is 
implemented by 



{Rt/m) VMD(SDL) 



{ml,mlj) 



(119) 



with 



VMD(SDL) 



To 



2 2 \ 



2/2 



l^cut 



1 + (72 + P2 

/ 2 2 2 



-Re^ 



™2 ,.2 

M, 



SM^ 



Mo 



SM^ 



+P2 



+ 



2/3 ' 



Mo 
SMy 

Mo 



Mo 



{mi,mM,M^,) 



2 2 2 \ S.A^ V / 2 2 71 /tI 

mi , rriM, - -r-;-{rni , m^,, M^„ 



1 + (73 + P3 

l^cut 



-Rc< 



5Ma, 2 2 2 \ ^Ma, 2 2 7i^2n 



A4o 



Mo 



+(73 



+P3 



l^cut 



5Al4. 2 2 _ 



2 N SMa^ 2 2 .^2n 



5M 
Mo 



A , 2 2 2 \ 

'mi,mM,^^cut) 



Mo 

^(m^mi„M» 



(120) 



where /2 and /s are flags which are either one or zero, determining whether or not the 
photon coupling to the vector resonance and the pion or to the axial resonance and the 
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pion are dominated by uj and p type resonances, respectively. The parameters and p2 
for the 10 type ones and (T3 and for the p type ones describe the relative contributions 
of higher radial excitations. 

From our results we can also calculate the correction 5Re/fj, to the ratio of the electronic 
and the muonic decay modes of the pion: 

^Re/ii = i^Re/fj.) CPM + {SRe/^i)vMD(PML) + {SRe/^i)HSD 

+ {.5Re/n)vMD(HSD) + iSRe/ti)sd (121) 

where {5Re/^)cPM^ ■ ■ ■ , {SRe/iJ.)sd are defined and calculated in a completely analogous 
way. 



9 Numerical Results 

We will now evaluate the formulae of the last section numerically and so get our final 
results. For the calculation of the standard loop integrals we have used the numerical 



programs FF and AA/FF 

Uncertainties of the final result come from two different sources: on the one hand from 
uncertainties in the hadronics and on the other hand from uncertainties in the matching 
of long and short distances. We estimate their sizes by varying the various parameters 
involved, ie. FJ(0), relative contributions of higher radial resonances and others for the 
hadronics and by varying pcut for the matching. We obtain a central value with an error 
estimate in the following way. We discuss the various parameters and find for each one 
of them a central value and some reasonable range over which we will vary them. Taking 
all parameters at their central values, we obtain the central value for the total radiative 
correction. We then vary the various parameters and determine which choices give the 
smallest and the largest radiative correction, respectively. Then by simultaneously taking 
those values which result in the largest (smallest) correction, we obtain an upper (lower) 
limit for the radiative correction. 

For the parameters of the hadronics, we obtain in this way three parameter sets (I)- 
(III), we are given is the tables [^, ^ and H, corresponding to the central values (I) and to 
the lower (II) and upper limits (III) on the total radiative correction. 

For the scale pcut the range pcut = 1 ... 2 GeV is reasonable and we will use pcut = 
1.5 GeV as an intermediate standard value. 

If all parameters have their standard values we find the total radiative correction 

Rr/^ = 0.1Q% (122) 

Because of all the ambiguities of the non-perturbative strong interactions in the 
0(1 GeV) regime, it is interesting to consider the result which is obtained if in the long 
distance part only the point meson is used, without vector meson dominance and with- 
out hadronic structure dependent loops. This point pion contribution is fixed by the low 
energy theorems of QCD and therefore free from any ambiguities. However, in this case 



the result depends very strongly on the value of the scale pcut- This is shown in Fig. [15 
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Table 1: Parameter sets for i?^/, 



Paraiiioter 


I 


II 


III 




Fa{0) = 




0.0116 


0.0100 


0.0132 




[MeV] = 




400 


600 


250 




(7 = 




0.136 


0.0584 


0.000 




P = 




-0.051 


0.0 


0.000 


widths 


in real radiation: 


Fv 


energy dependent 


energy dependent 


fixed 


widths 


in real radiation: 


Fa 


energy dependent 


fixed 


energy dependent 








-0.1 


-0.1 







Pi = 




-0.04 


-0.04 







/2 = /s = 




1 





1 




cr2 — 0-3 — 









-0.1 




P2 = P3 = 









-0.04 



Table 2: Parameter sets for Rr/x 



Parameter 


I 


II 


III 


F^iO)- 


0.0525 


0.0425 


0.0625 


Tk, [MeV] = 


90 


110 


70 


ax = 


0.000 


0.136 


0.000 


Pk = 


0.000 


-0.051 


0.000 


widths in real radiation 


fixed 


fixed 


fixed 


(Ti = 





-0.1 





Pi = 





-0.04 





f2 — h — 


1 





1 


C^2 = 0"3 = 







-0.1 


P2 = Pa = 







-0.04 
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Table 3: Parameter sets for Re/^ 



jT ai aiiietei 


T 
i 


TT 
11 


TTT 
111 


^A(Uj - 


U.Uiio 


U.UiUU 


n ni Qo 


"p rT\/r,^"\7'i 
i ai [MeVJ = 


4UU 


^^nn 
oOO 


o txn 

zoO 


(T = 


U. iou 


n nc;8/i 


n nnn 
u.uuu 


P = 


— U.UOi 


n n 

U.U 


n nnn 
u.uuu 


VV1L4.L110 111 IC-CXl 1 CHJ-lCl) tlv^ll. -1 y/ 


pTip-rcTAr nPFtprinpnl" 


prip-rcrv HpTtPTinPnt 
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where we show the correction ^-Rt/tt in variation with the scale ficut both for our complete 
calculation, using standard parameters, and for the calculation where in the long distance 
part only the point meson is taken into account. In the last case the dependence on /Xc^t 
is very strong. As has been discussed before, the pure point meson is reliable only for 
Pcut ^% t)ut the short distance correction only for ficut ^ (1 • • • 2) GeV. So here a very 
large range for ficut has to be considered. For a small ficut the correction 6Rr/n could 
even become negative, and at ficut = 2 GeV the correction is about SRr/n = 0.7% and 
still rising strongly. 

This large dependence on ficut results from the incomplete treatment of the long 
distance part and hints to the necessicity of improving the model. The inclusion of VMD 
and of the structure dependent loops decreases the dependence on ficut very much, as 
can be seen from the other curve in Fig. O. Above ficut = 2 GeV the curve becomes 



almost completely fiat, and the variation in the (relevant) range ficut = (1 • • • 2) GeV is 
smaller than ±0.05%. In Tab. ^ we show the different contributions which add 

up to the total radiative correction, using ficut = 1.5 GeV. The individual contributions 
involve moderately large logarithms such as ln{mr/m.^) or ln{mp/mr). However they 
have opposite signs, such that most of the corrections cancel and only a very small total 
radiative correction of SRr/n = +0.16% is left. 

In Fig. |TB| we show the long and short distance corrections in variation with ficut- 
We display separately the correction due to the effective point meson with vector meson 
dominance coupling in the photon coupling 

{SRt/m) cpm + {SRt/m) vmd(pml) (123) 
the hadronic structure dependent correction with vector meson dominance 

{SRt/m)hSD + {SRr/M)vMD(HSD) (124) 

and the short distance correction 

{5Rr/M)sd (125) 
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Tabic 4: The different contributions adding up to tfie total radiative correction SRt/tt for 
l^cut = 1-5 GeV 



effective point pion 




-|-i.UO /o 


cutting off tlie point pion loops 

&t l^cut 


i 


-0.21% 


vector meson dominance of the 
pion electromagnetic form factor 




-0.38% 


Structure dependent radiation 
and SD-IB interference (real 
photon emission) 


^sd+int{t ■nvrl) 

Vq{t TTVr) 
ro(7r ^.U^) 


+0.05% 




5^ SDL, 2 2 n\ 

(m^,m^,0) 

i 




hadronic structure dependent 
loops, cut off at ^cut 


J- 

(to TO^,0) 

to 

i 


-0.24% 


vector meson dominance of pho- 
ton couplings in the hadronic 
structure dependent loops 


{SRr/n)vMD(SDL) 


+0.13% 


short distance contribution for 


{SRr/n)sd 


-0.25% 


sum 




+0.16% 
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It can be seen clearly that the long and the short distance corrections vary in the opposite 
way with fj,cut, such that the dependence of the total radiative correction on fi^ut is 
significantly smaller than that of the individual corrections separately. This is a sensible 
result. If our long and short distance corrections would exactly describe the real world 
within an overlapping region of fj,cut, the sum of the two would be independent of ficut 
within this overlap region. The small remaining dependence of our final result on ficut 
within the range ficut = 1-0 . . . 2.0 GeV indicates that our model is not unreasonable. 

T7\ we show 6Rr/n in variation with ficut for three different choices for the 



In Fig. 



parameters of the resonance physics, viz. for the standard set (I) and for the set (II) 
which gives the smallest correction and for the set (HI), which gives the largest correction 
(see Tab. |I|). All the three curves have been obtained using real vector meson masses 
(narrow width approximation) in the virtual corrections, but we have compared with 
the results obtained with complex vector meson masses, and in all cases the difference is 
extremely small and completely negligible. While the two curves for (I) and (II) are close 
together, the curve for (III) lies significantly higher. The large difference is due to the 
question whether or not the photon in the hadronic structure dependent loops couples 
to the mesons via vector meson dominance, all other parameter variations have a much 
smaller influence. 

As our flnal result on 6Rr/n we get from the parameter sets (I) — (III) and from 
varying ficut = 1.0 ... 2.0 GeV: 



6Rr/n 



0.16" 



-0.09 
-0.14 



% 



(126) 



In Fig. |1^ we show 5Rr/K in variation with ficut for the three parameter sets (I)-(III) 
which are defined in Tab. |^. Our final result for the correction to Rt/k is 



SRr/K = (0.90l[j-]^) 



(127) 



where the central value is for the parameter set (I) and with ficut = 1.5 GeV, and the 
lower and upper limits are from the sets (II) and (III) and ficut = 1.0 GeV, 2.0 GeV, 
respectively. 

For the normalized branching ratios this results in 



expt 



theo 



_ n 61 99 + fr0007 +0-0005 _n«ioq+0.0009 
BR~ ~ 0.6129 ± 0.0007 _Q QQQg -0.6129_Q QQ^;^ 



BR 



expt 



theo 



= 0.0406 ± 0.0002 ioooo5 = 0.0406 ± 0.0002 



expt 



theo 



= 0.6535 ±0:^i0;0005 = 0.6535+0-0009 



BR, 



-0.0012 



'128) 



The first errors given (called "experimental" ) are due to the uncertainties in the lifetimes 
of the mesons and in the tau mass m,-, and the second errors (called "theoretical") are 
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due to the uncertainties in the radiative correction. This final result deviates from the 



experimental result p6| given in the introduction 

/BR^^ ^(11.99 ±0.25)% 

V BRe /e^p (17.76 ±0.15)% ^ ' 

by 1.4 standard deviations. The agreement between theory and experiment is not signif- 
icantly enhanced by the inclusion of the 0{a) corrections to the decay rate for r ttUt-. 
Thus if the standard model is correct, either the current experimental number for the 
branching ratio BR^+x is slightly too large or the one for BRg slightly too small, or both. 

Some comments on the reliability of our matching procedure are in order. First it is 
important to note that indeed the dependence of the long and short distance corrections 
separately on ^cut is considerably larger than that of their sum, as has been discussed 
above in connection with Fig. |16[ Second we have performed the matching in a certain 
way be splitting the photon propagator into a long and a short distance part according 
to Eqn. (|35D. This corresponds to a soft transition from long to short distances. Another 
way to perform the matching would be a sharp transition from long to short distances 
at Hcuti by integrating the long distances from k\ = 0. . ■ fil^f and the short distances 
from /c|; = /iciit • • -rn^z- principle the results obtained with this method might differ 
somewhat from our results. However, the uncertainty of our final result for 6Rr/n is 
dominanted by the hadronic uncertainties and not by the matching uncertainties. Varying 
the parameters of the hadronics and the matching scale ficut separately, we obtain 

/ matching hadronicsX 



SRt/tt 



n +0.02 +0.07 
^■^^ -0.05 -0.13 



(130) 



V / 

which clearly displays the dominance of the hadronic uncertainties. Therefore we think 
that the precise procedure for performing the matching is not essential. Still it might be 
an interesting task to repeat the calculation with a matching based on a sharp transition 

fJ-cut- 

Now we will use our results to predict for the pion decay the ratio Re/^ of the electronic 
and muonic modes. In Fig. |l^ we show 6Re/^ in variation with ^^ut for three parameter 
sets (I) — (III) defined in Tab. We find for the total radiative correction 

5R^f^ = (-3.79 ± 0.01) % (131) 

But note that this error, which is due to the uncertainties in ^cut and in the hadronic 
resonance physics, is smaller than the one which is to be expected because of the neglect 
higher order corrections of O(a^), as we will explain below. 
The value in Eqn. (|131|) results in 

Re/, = <^l + 5i?e/,; 
i2 




1 + 6R, 



1.2835 ■ 10-^(l + (-3.79 ± 0.01) % 
(1.2349 + 0.0001) ■ 10"^ (132) 
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where the error given includes only the uncertainties resulting from the hadronic structure 
dependent corrections and the matching scale /icut and misses the comparable uncertainty 
resulting from the neglect of the 0{a) corrections. 

This result is to be compared with the recent result given by Marciano and Sirlin [H: 



= (1.2352 ± 0.0005) ■ 10"^ (133) 



which is equivalent to 



SR^/^ = (-3.76 ± 0.04)% (134) 

Their calculation is rather similar to ours but differs in some details. In the short dis- 
tance corrections they use the leading logarithms only, supplemented by leading QCD 
corrections and furthermore by summing up the leading logarithms to all orders in a 
via the renormalization group. But since the short distance corrections cancel almost 
completely in Re/fj., these differences are not important. For the long distance corrections 
they also consider both the effective point pion field and hadronic structure dependent 
corrections. However, for the latter they do not use a specific model but rather use a 
theorem given by Terent'ev P7| on the leading lepton mass dependent hadronic structure 

effects of the order ^In — ^. The possible size of the remaining hadronic structure 

vr iTLp mf 

dependent corrections, which they do not calculate, represents the main source of the 
uncertainty in their Re/^ as quoted above. The authors also consider the effect of cor- 
rections of higher order in a. Since the correction 3(a/7r) ln(m^/me) ~ —3.7%, which 
results from the lepton mass singularities, dominates the 0{a) correction, one expects 
its higher order counterparts to similarly dominate their respective orders. Summing all 
such logarithms via the renormalization group gives the enhancement factor 



1 _ |£i In !IM 

in me 



1-^ln^ 



9/2 

— = 1.00055 (135) 



which multiplies Re/fj.- This implies that the O(a^) correction to R^/^ is of the order of 
0.05%. 

If we multiply our result with this enhancement factor (see Eqn. ( |135| )), we get 

= (1.2356 ± 0.0001) ■ 10"^ (136) 

or 

5R^/^ = (-3.74 ± 0.01) % (137) 
which agrees with the Marciano-Sirlin result within their error bars. 



10 Summary and Conclusions 

We have calculated the radiative corrections to the decays r Mvj. and to M lh>i 
[M = TT OT K, I = e or fi). The central issue of this paper was the treatment of the strong 
interaction. The amplitude of the radiative decays with the emission of a real photon 
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can be divided into the amplitude for internal bremsstrahlung (IB) and the amplitude 
for structure dependent radiation (SD). In the virtual corrections one has to integrate 
over the momentum k of the virtual photon, and therefore all the three energy regimes of 
the strong interaction have to be taken into account. For small k'^ we have contributions 
from the long distance correction, which consists of the point meson contribution and of 
the hadronic structure dependent part. The short distance corrections, which we have 
calculated using the parton model, contribute in the large fc^ region. 

Our final result for the radiative correction 6Rr/n to the ratio r(r 7rz/T-(7))/r(7r 

6Rr/^ = (o.ielg-OQ) % 

and for the ratio r(r — > Kh'T-{'y))/T{K — > yUz/^(7)) we obtain 

5R^/K = (0.90+0;^) % 

Note that these numbers are calculated by summing up virtual, soft and hard photonic 
corrections. 

We can translate the radiative corrections into predictions for the branching ratios 
BR(r ^ 7rz/^(7) = (11.10 ± 0.02)% x ^ 



BR{T-^Kur{l) = (0.737 ± 0.005)% X 

BR{t ^ huri-f) = (11.84 ± 0.02)% X f — (138) 



295.7 fs 

295.7 fs 
295.7 fs 

where h denotes the inclusive sum of pions and kaons. 

For the ratio Re/^ of the electronic and muonic decay modes of the pion, we obtain a 
radiative correction of 

SR^/^ = (-3.74 ± 0.01) % (139) 

resulting in 

Re/f, = (1.2356 ± 0.0001) ■ 10"^ (140) 
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Appendix 

Except for the tau mass we use the standard particle data of [|12] for masses and widths. 
In the case of the tau mass, we use 

= (1777.1 ± 0.5) MeV (141) 
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In this paper we use the notation of Bjorken and Drell |]28|, especially 

eoi23 = +1 7' = ^7°7SV (142) 
and the e is the charge of the electron, e < 0. The standard loop integrals are defined by 



167r2 

4_£) / d k 

/i 



[Bq I B^ I 5^^] (pi , mo , mi 



(27r)« P - ml 
4—/) /" dPk [1 I A;^ I k^ky] 



167r2 



(27r)^(P-m§)p + pi)2-mf] 
d k [1 I I k^ki,] 



(27r)^ - mg) [(A; + pi)2 - mf] [(A; + ps)' - mi] 



(143) 

The scalar functions depend on invariant combinations of the momenta only, 

5o(pi,mo,mi) = 5o(pi,mo,mi) 
Co{pi,p2,mo,mi,m2) = Co(Pi, (pi - P2)^ "^o, "^1, mg) (144) 

and the vector and tensor integrals are decomposed covariantly in the form 

B^{pi,mo,mi) = BipI 

5^'^(pi,mo,mi) = Boog^" + BupIpi 

C''{pi,p2,mo, nil, 1712) = C1PI + C2P2 

C^'^(pi,P2,mo,mi,m2) = Coo^?'^^ + CupW + C22P^P2 + ^^12 (rfPs + (145) 



where 



= B,iplmo,mi) 

E,, = i?,,(p?,mo,mO 

Ci = Ci{pl,{pi-p2Y,pl,mo,mi,m2) 

Cij = Cij{pl,{pi - p2f,pl,mo,mi,m2) (146) 
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Figure 1: The three different classes of radiative corrections 



Figure 2: Short distance corrections 



Figure 3: Long distance corrections 



Figure 4: Point pion diagrams 



Figure 5: Short distance diagrams 



Figure 6: Decay r vri/^ via an intermediate quark-antiquark state. The bubble on the 
left hand side is to be replaced by the short distance diagrams and 6Ai, cf. Fig. |^ 



Figure 7: The function rT-/j,{u) for fi^ut = 1 GeV (dashed), ficut = 1-5 GeV (dotted) and 
ficut = 2 GeV (solid) 
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Figure 8: The function r^jKiu) for Hcut = 1 GeV (dashed), Hcut = 1-5 GeV (dotted) and 
Hcut = 2 GeV (sohd) 



Figure 9: The function re/^{u) for ^cut = 1 GeV (dashed), ^cut = 1-5 GeV (dotted) and 
^icut = 2 GeV (sohd) 



Figure 10: Short distance correction {SRr/Tr)sd- complete result according to Eqn. 
(dotted) and estimate based on the leading logarithms according to Eqn. (03) (solid) 



Figure 11: Diagrams for hadronic structure dependent corrections 



Figure 12: Vector meson dominance of coupling of the photon to the pion 



Figure 13: Vector meson dominance of the photon coupling in the hadronic structure 
dependent diagrams 



Figure 14: Vector meson dominance in the seagull coupling 
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Figure 15: Radiative correction to Rt/-k'- Complete prediction using the standard param- 
eter set (solid) and prediction for short distance plus point pion only, ie. {SRr/n)cPM + 
{SRT/-K)sd (dashed) 



Figure 16: The different contributions which add up to the total correction: The point 
pion correction, supplemented with vector meson dominance in the photon coupling 
(dashed), the hadronic structure dependent correction with vector meson dominance in 
the photon couphng (dash-dotted), short distance correction (dotted) and their sum, the 
total correction (solid) 
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Figure 17: The total radiative correction SRr/i^ for different choices for the parameters of 
the structure dependent correction: Standard choice (I) (sohd), choice (II) (dashed) and 
choice (III) (dotted) 



Figure 18: The total radiative correction Rt/k for different choices for the parameters in 
the structure dependent correction: Standard choice (I) (solid), choice (II) (dashed) and 
choice (III) (dotted) 



Figure 19: The total radiative correction Re/fx for different choices for the parameters in 
the structure dependent correction: Standard choice (I) (solid), choice (II) (dashed) and 
choice (III) (dotted) 
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